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Abstract 

We study theoretically the chirality of a generic rigid object's sedimentation in a fluid under 
gravity in the low Reynolds number regime. We represent the object as a collection of small Stokes 
spheres or stokeslets, and the gravitational force as a constant point force applied at an arbitrary 
point of the object. For a generic configuration of stokeslets and forcing point, the motion takes 
a simple form in the nearly free draining limit where the stokeslet radius is arbitrarily small. 
In this case, the internal hydrodynamic interactions between stokeslets are weak, and the object 
follows a helical path while rotating at a constant angular velocity ui about a fixed axis. This u 
is independent of initial orientation, and thus constitutes a chiral response for the object. Even 
though there can be no such chiral response in the absence of hydrodynamic interactions between 
the stokeslets, the angular velocity obtains a fixed, nonzero limit as the stokeslet radius approaches 
zero. We characterize empirically how uj depends on the placement of the stokeslets, concentrating 
on three-stokeslet objects with the external force applied far from the stokeslets. Objects with the 
largest u are aligned along the forcing direction. In this case, the limiting uo varies as the inverse 
square of the minimum distance between stokeslets. We illustrate the prevalence of this robust 
chiral motion with experiments on small macroscopic objects of arbitrary shape. 

PACS numbers: 47.57.ef, 47.57. J-, 87.16.Ka, 47.63.M- 
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I. INTRODUCTION 



It is not unusual to see objects falling through water or air twisting as they sink. For 
example, a propeller-like maple seed will twirl as it falls from the tree. A consistent preference 
for twisting in a particular direction would constitute a chiral response of the object. Such 
a response must reflect some chirality in its shape, and the magnitude and nature of the 
twisting is evidently a consequence of well-known hydrodynamic laws. However, there is 
little fundamental understanding of what features of the shape control the magnitude of a 
chiral response. 

In the past decade there has been a revival of interest in the tumbling motion exhibited 
by extended objects as they fall through air [1, 2]. These complex motions are of a different 
nature than what we study here. The objects under consideration have no intrinsic chi- 
rality, and interesting motions depend instead on significant Reynolds numbers, where the 
advection of momentum through the fluid is important. 

Aside from these, a few studies have examined the low Reynolds number sedimentation 
of different bodies. For a specific propeller-like design, Makino and Doi [3] showed that an 
ensemble of identical particles with different initial orientations will bunch together into a 
cylindrical shape oriented along the direction of the sedimenting force, whereas a similar 
group of achiral ellipsoids will drift apart. They have also made some headway in classifying 
the range of allowable motions for objects depending on whether or not they are skew or if 
there is an applied torque [4]. Gonzalez, Graf, and Maddocks [5] have further explained some 
properties of the possible motions. We hope to improve on the parts of this understanding 
related to chiral objects. 

Understanding the connection between shape and chiral motion would allow chiral sedi- 
mentation to be used as a characterization tool for objects of a supramolecular scale, such 
as colloidal particles and cells. Detecting the rotation of sedimenting bodies would give 
information not obtainable from other simple probes such as dynamic light scattering and 
intrinsic viscosity. These conventional measures sense only the hydrodynamic size of the ob- 
jects, whereas rotation speed can sense the distinctive feature of a permanent chiral shape. 
Many biological structures have a strong chirality that is unrelated to propulsion. Examples 
include protein-DNA complexes [6, Chapter 10] and fibrils such as Actin [6, Chapter 22], 
which are made of repeating subunits. Such objects must rotate as they sediment, and a 
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knowledge of the connection between their shape and their rotation would be valuable. 

We will show that chiral motions are natural to characterize when the hydrodynamic 
interactions between parts of the body are small. Thus, much of our study will be aimed at 
objects with this property, which we will term "nearly free draining." Physical realizations 
of such objects can include thin, rod-like objects such as microtubules [6, Chapter 23], 
bacterial flagella [6, Chapter 23], or sickle cells [7, Chapter 19]. As a concrete example, we 
can consider the propeller shape of Makino and Doi [3] , shown in Figure 7 below. For such 
an object of length about 10 microns in water, with a density of about 1 g/cm 3 , we predict 
rotational velocities on the order of 10 Hz. This should be noticeable, even when compared 
to the rotational diffusion coefficient, which for an object of this scale is only of order 10~ 4 
Hz. Smaller objects on the scale of a micron or less will also have a noticeable effect if they 
sediment under a slightly larger force, as in a centrifuge. 

In Section II, we discuss the equations of motion for our objects, and show how any 
inherent chirality must be encoded and expressed. In Section III we introduce the "tumble 
zone," a region in parameter space which determines whether or not a sedimenting object 
can exhibit ongoing tumbling behavior, and put a bound on its size. Following that, in 
Section IV we review a stokeslet formalism for modeling rigid bodies, and show how to use 
this to calculate the internal hydrodynamic interactions needed in the equations of motion. 
In Section V, we find these interaction effects in the nearly free draining limit, where the 
interaction strength becomes small. In this limit, we find that the tumble zone becomes 
arbitrarily small, and that almost all objects will exhibit chiral sedimentation. Once this 
is established, in Section VI we show how the chiral response behaves in certain limiting 
cases. Using a simple three-stokeslet body, we empirically examine how different aspects of 
shape affect our measurement of chirality. In Section VII we show the results of numerical 
simulations. We check these numerical results against the analytic ones found in Section V, 
and compare the typical motions of a random chiral body with both a more symmetric 
propeller shape and an achiral ellipsoid. Finally, we report the results of a simple experiment 
done on small macroscopic objects of arbitrary shape. 

Throughout the next several sections, we refer to many different types of objects. To 
distinguish them, we use the following conventions: 3-vectors and unit 3-vectors will be 
denoted with arrows (e.g. v) and hats (-0), respectively. The 3x3 matrices that operate on 
them will use a blackboard bold font (M). 6- vectors will be in italics with vector signs (V), 
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and the 6x6 matrices will be underlined (M). Large vectors composed of 3- vectors for each 
stokeslet will be bolded with vector signs (v), and the matrices that interact with them will 
be bolded with underlines (M). 

II. THE PROPULSION MATRIX 

In order to analyze the behavior of our sedimenting body, we take advantage of the 
fact that at low Reynolds numbers, the force and torque on a body are proportional to its 
velocity and angular velocity. Following Purcell [8], we collectively refer to these constants 
of proportionality as the propulsion matrix P. That is, we define extended force and velocity 
vectors T = (F, r) T and V = (V, uj) t , and write 



As a consequence of the Onsager relation, and the requirement that the dissipated energy 
be positive, this propulsion matrix must be both symmetric and positive-definite [9], so it 
can be written in block form as 



where IK and fl are symmetric 3x3 matrices which are also positive-definite. 

The propulsion matrix contains all of the information necessary to describe the dynamics 
of the object. Once it is known, an analysis of the motion can be carried out without 
reference to the specifics of an object's shape. 

In order to specify a torque, P must be computed about a specific point. Moving this 
point will change both C and fl, though IK will remain the same. Happel and Brenner [9] 
show how each of these individually transform under a change of coordinates. We arrive at 
equivalent results in a slightly different form. To begin, let £ and E' represent two different 
inertial frames used to describe variables. In the following, primed variables will denote 
quantities viewed in the E' basis, and unprimed ones will be those living in the E basis. We 

— » — # — # — * 

then have propulsion equations for each of the frames: T = JPV and T' = JP V . 

It is easy to transform between coordinate systems that differ only by a rotation: if R is 
the rotation matrix that will take one set of axes to the other, then each subblock X of P 
changes as X — > RXR -1 . Next we consider frames E and £' which differ only by location of 



f = PV 



(1) 




(2) 
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— * —* 

the origin, and let R be the vector to E"s origin. We now consider the effects of a force F 
and torque f applied at the origin of E. The body will feel the same net force and torque, 

— * — * 

and thus respond with the same motion, that it will if we pull at R with force F and supply 
a torque of f + (-R) x F. That is, 

f = f+ \ ® _ I = (1 + B)^, (3) 




-R) x F 

where the matrix B is defined in block form by 

\-lRxj 

Here we are using the notation that for any vector X, [Xx] is the antisymmetric 3x3 
matrix which satisfies [Ix]t) = X x v for all vectors v. 

There is some extended velocity vector associated with this given force and torque, but it 
will be represented differently in E and E'. The angular velocity must be the same in both 

— * — * — * 

systems, but a different linear velocity needs to be used. Using V' + Rxu = V + Oxuj,we 
can conclude 

(-R x uj \ t _ -, 
= (1- I T )V. (5) 
/ 

We can now combine these two expressions to get a relationship between JP and £': 

P = £'V' 
(1 + B)/ = P'(1-B T )V 

PV = (1 + B)"V(1- B T )V. 

— * 

Since V is an arbitrary velocity, we can just write 

£ = (I - B)P'(1 - B T ) (6) 

where we have used the fact that (1 + B) -1 = (1 — B). 

There is a unique point, termed the "center of reaction," [9] about which the sub matrix 
C is symmetric. For many objects with a high degree of symmetry, this often coincides 
with the centers of mass and buoyancy, but for a general case, these different points are not 
related. 
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For the sedimentation processes that concern us, it is often convenient to deal with the 
inverse of the propulsion matrix, known as the mobility matrix ]tf . We will write it in block 
form as 



The matrix A, which we will call the alacrity matrix, gives the velocity response to an 
applied force. Our screw matrix S gives the angular velocity caused by supplying a torque, 
and the twist matrix T shows the coupling between angular velocity and force. 

Since P, IK, and fl are symmetric and positive definite, M, A, and S must be as well. 
There is also a unique choice of origin for which T is symmetric, but it is in general different 
from the center of reaction. We will call this point the center of twisting. Furthermore, by 
inverting the transformation law for P, we can find how M changes if the origin is moved: 



From this, one can see that the screw matrix remains invariant (S = S'), and that the twist 
matrix changes simply as 



Conveniently, for sedimentation processes the twist matrix alone captures both the chiral 
information and the dynamics of interest. Indeed, if we want our sedimenting object to 
show a preferred chirality, jM must not be invariant under inversions about the origin. A is 
necessarily invariant under this inversion, since both force and velocity transform as vectors. 
The same is true for the screw matrix, since both torque and angular velocity transform 
as pseudovectors. However, the twist matrix will reverse sign. Thus any chirality in the 
object must manifest itself through this twist matrix. As a simple example, we see that if 
the center of twisting is at the origin, then an object can only be chiral if the eigenvalues of 
T are not symmetric about 0. 

The physical manifestation of chirality we are concerned with is the rotation of our object: 
at any time t, u(t) = J(t)F + S(t)r. However, since sedimentation involves forces acting 
on the centers of both mass and buoyancy, with no supplied torques on either, it is possible 
to choose as our origin a point of zero torque on the object. In this case, we just have 

— * 

uj(t) = J(t)F, which allows us to restrict our attention to the twist matrix. 




(7) 




(8) 



T = T' — Sj-Rx]. 



(9) 
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The twist matrix scales in a simple way with the overall size of the sedimenting object [9]. 

— * 

For a given object, the force F needed to produce a given uj is proportional to the viscosity, rj. 
Thus T is inversely proportional to rj: J = t]^ 1 !, with T independent of viscosity Evidently 
T has dimensions of viscosity / (force x time), or (length) ~ 2 . With a fixed force, the rotation 
rate for an object enlarged by a factor a will thus be reduced by a factor of a 2 . Analogous 
reasoning shows that the velocity V is reduced by a factor a 3 . The object's translation for 
a given increment of rotation thus varies linearly with a, and enlarging the object simply 
enlarges the path of its sedimenting motion by the same factor. 

Given the twist matrix at time t, it can be found some small At later by rotating T(t) 
by the angle uj{t)At. Then 

J(t + At) = (l- At[w(i)x])T(t)(l + Af[w(i)x]). 

Eliminating terms of order At 2 gives 

J(t + At) = T(t) + Af[T(t), Mt)xJ] 

Taking At -> yields 

t=[T,[TFx]]. (10) 

The evolution of the other blocks of the mobility matrix evolve in a similar fashion: 
A = [A, [TFxJ], and likewise for S. 

This formalism, with fixed axes in the lab frame and a dynamical T, is equivalent to the 
Euler equation formalism used by Gonzalez et. al. [5], which treats the body axes as fixed, 
and considers a dynamic force vector. We denote quantities in this body frame of reference 
using double-prime marks, ". At each instant the body frame rotates relative to the space 
frame with angular velocity Cu, as noted above. Thus the space frame rotates with respect 
to the body frame at angular velocity —uj, and F" = —uj x F" . This uj, common to both 

— * — * 

frames, can be expressed equally in the body or space frame: uj = JF = J"F". 

Of particular interest are stationary states, in which the essential part of the motion is 
constant in time. In the body frame, a stationary state is one in which F" = 0. Since 

F" = -wx F" = —J"F" x F", (11) 

— # 

there is a stationary state if and only if F" is an eigenvector of T", with an eigenvalue that 

— * — * 

we denote as A. Since F" is constant in time, u = J"F" must be as well. 
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Because the twist matrix is 3 x 3, it has either one or three real eigenvalues. In the 
case of a single real eigenvalue, the analysis above implies two fixed-point forces in opposite 
directions. The sign of the eigenvalue gives the chirality: a positive eigenvalue means that 
with the usual right-handed definition of angular velocity, the object twists as it descends 
in the direction of a right-handed screw. The chirality of the two fixed points is thus the 
same. However, the stability is not. The stability of the fixed point direction Fq can be 
determined by considering the quantity F" ■ Fq. Its derivative F" ■ Fq determines whether 
F" moves toward or away from the fixed point with time. One may readily show [5] that for 
a given Fq, the sign of this derivative is fixed for all F" ^ ±Fq. If this were not the case, 
then there would be some F" for which F" • Fq — 0. To see that this is impossible, note that 
it either requires F" _L Fq or F" = 0. Consider first the case where F" _L Fq. Equation 11 
tells us that F" _L F" and F" _L J"F". Since Fq is the only eigenvector of T", this means 
that F" has no component in any of the three independent directions Fq, F", and J"F". 
This leaves us with the option that F" = 0. However, this means that F" is a fixed point, 
which contradicts the assumption that Fq is the only eigenvector of T". Thus F" ■ Fq ^ 
for all F" 7^ ±Fq, meaning that F" ■ Fq has the same sign for all such F" . If the sign is 
positive, then all F" move toward the Fq axis and Fq is then a globally stable fixed point. 
Evidently the opposite fixed point at —Fq is globally unstable. 

With three real eigenvectors, F" ■ Fq can vanish at points besides ±Fq, so the global 
stability argument above is no longer valid. The simple chiral signature of the object is no 
longer present, and the motion becomes more complicated and depends on initial conditions 
[5]. Happily, this case can be excluded for a large class of objects, as we show below. 

III. THE TUMBLE ZONE 

Given a fixed shape for an object, we can choose the center of twisting as our origin. At 
this point the twist matrix is symmetric, meaning T must have three real eigenvalues. Next, 
keeping the object's shape fixed, we can explore the locus of points to which we can move 
the forcing point while still keeping all three eigenvalues real. Since, as shown above, global 
stability is not present at these forcing points, we call the region they form the "tumble 
zone" for that particular shape. We will show here that the volume of this tumble zone is 
always finite. 
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With any choice of origin, the screw matrix S is always symmetric with positive eigen- 
values, as discussed in Section II. We may then work in the basis where 

S = diag(si,s 2 ,s 3 ). 

In this basis, we will move the forcing point to R p . From Equation 9, this will give us a 
new twist matrix T = V — S[-R p x], where TP is the twist matrix computed about the center 

— * 

of twisting. We will show that if we choose R p to be sufficiently large, then the new twist 
matrix about this origin must have only one real eigenvalue. 

We can compute the discriminant A of the characteristic polynomial of our new T. If 
the discriminant of a cubic equation is positive, then there is one real root and two complex 
conjugate ones. In this case, our twist matrix will have only one real eigenvalue. The 
discriminant is 

A = 27Det 2 (T) - 4Det(T)Tr 3 (T) + 9Det(T)Tr(T)(Tr 2 (T) - Tr(T 2 )) 

-iTr 2 (T)(Tr 2 (T) - Tr(T 2 )) 2 + -(Tr 2 (J) - Tr(T 2 )) 3 , (12) 
4 2 

which is homogeneous of degree 6 in T. 

The discriminant A is a sixth degree polynomial in R p , so if the coefficient of the R® 
term is positive, we can be assured of getting A > for any R p bigger than the largest 
root of this polynomial. Since A is homogeneous, there can be no powers of ¥' in the R p 
term. This leading term can thus be found from Equation 12 by replacing T with Sp? p xJ. 
Since S is symmetric and [i? p x] is antisymmetric with a zero eigenvalue, Det(S[i? p x]) and 
Tr(S[i2pX]) both vanish. Accordingly, the only term in Equation 12 that can contribute in 



order R® is the last one: 



A = -i 
2 



3 



Tr[(S[i? p xJ) 2 ] + 0(R 5 p ) (13) 



In terms of the eigenvalues Sj and the coordinates R p i, R P 2, and R p3 , this trace has the form 

Tr[(S[i? p xJ) 2 ] = -2(R 2 pl s 2 s 3 + R 2 p2 s lS3 + R 2 p3 s lS2 ). 
Since the Sj are all positive, if we define s m = minjsj} then we can write 

A>2s e m Rl + 0(R 5 p ), (14) 
whose leading term has a positive coefficient. 



Thus outside a sphere of sufficient radius R p the discriminant is positive, there is a 
single real eigenvalue, and the motion converges to the globally stable motion discussed in 
Section II. 

IV. THE STOKESLET REPRESENTATION 

The propulsion matrix for a body can sometimes be found analytically, and there are 
several known results for objects with various symmetries [9]. However, it can be more 
difficult to find when such symmetries are not present. We use the approach of Kirkwood and 
Riseman, as described by Meakin and Deutch[10], in which a sedimenting body is represented 
as a rigid collection of small beads known as stokeslets. Each stokeslet corresponds to a 
point source of drag, which exerts a force proportional to its velocity: F = —^yv, with 
drag coefficient 7 = Qnrjp proportional to the fluid viscosity 77 and effective radius p of the 
stokeslet. 

By arranging the stokeslets appropriately, the flow field from most objects can be recre- 
ated [11]. Thus, they form a simple way of modeling arbitrary bodies. This approach is used, 
for example, to model flagellar propulsion [12]. Carrasco and de la Torre [11] investigate the 
effectiveness of different strategies for placing the stokeslets. 

To create a propulsion matrix from a collection of stokeslets, one must take into account 
the change in fluid velocity past each stokeslet caused by the presence of the others. If one 
does not include these hydrodynamic interactions, then there can be no chiral effects in the 
sedimentation: the object will sink straight down, so all drag forces will be vertical in order 
to oppose it, and thus there can be no torque about the vertical axis. However, if we do 
include these interactions, the velocity at each stokeslet may be perturbed from the vertical, 
possibly causing a torque about that axis. This can make the object demonstrate chirality 
by spinning. 

The tool we use is the Oseen equation, which gives the change in fluid velocity caused 
by one of these stokeslets. In the frame of a body with n stokeslets, let v be the 3n 
component vector containing the velocity of the fluid at the locations of each stokeslet, 
taking hydrodynamic interactions into account: v = (V 1 ,^ 2 , . . . ,v n ) T . We define F to be 
the 3n component force vector acting on the stokeslets, and v e the external (undisturbed) 
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velocity of the fluid at the location of each stokeslet. Then we can write the Oseen equation, 

v = v e + LF, 



where L is the Oseen tensor [13]. If we denote particle number by Greek letters, and cartesian 
coordinate by Roman letters, then for a ^ (3 we can write 



j a/3 _ 



1 



(15) 



with rf the ith coordinate of particle a, and r a ^ the distance between particles a and (3. 
For a = (3 we should have 0, since an individual stokeslet cannot affect itself. 

Let U be the 3n x 6 matrix which relates the 3n dimensional v e and the extended velocity 
vector V = (V,uj) t : v e = UV. Since the velocity of the fluid past each stokeslet is the 

— * 

opposite of the velocity at which the object is moving through the fluid, v° = —V — uxr 01 . 
Thus we can see that U = (U 1 , . . . U n ) T , with 



(-1, [Px] 



/ 



V 



-10 
0-10 



z y 

r" — r\ 




This U matrix also has the property that T = U F. 

If we define the 3n x 3n matrix T =diag(7 1 , 7 1 , 7 1 , 7 2 , 7 2 , 7 2 , 
and we can rewrite the Oseen equation as 

v = v e + LTv 



,7 n ,7 n ,7 n ), thenF = £v, 



v 

u T rv 

f 



fl - LT) 



-i--t 



Lfrd-LTJ-Ve 
u T r(i - LT) _1 UV 



But this is just our definition of the propulsion matrix: 



p = u T r(i - Lr) _1 u 



(16) 



This result, which is a straightforward extension of the Kirkwood Riseman method ex- 
plained in Reference [11], shows that one can calculate P from a matrix inversion. 
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V. THE NEARLY FREE-DRAINING LIMIT 



A. The propulsion matrix 

When using the stokeslet model, we have an obvious mechanism by which we can model 
nearly free draining bodies: we simply take the stokeslet size (and thus the drag coefficient 
7) close to zero. If we assume from now on that each stokeslet has the same effective radius, 
we can obtain perturbative expansions in this common 7, and write 

IK = [K07 + tK l7 2 + . . . 
fl = n 7 + !Hi7 2 + • • • 
C = C 7 + C17 2 + • • • 

To first order in 7, there are no hydrodynamic interactions, so we just have 

— # — # 

v = v e 

when the body is not rotating. In this case, the total force on the object is just the sum of 
the individual forces acting on each stokeslet: F = Yla=i(~ 7^ a ), and v a = —V is the same 
for all a. But F = (IK07) V, so we get a IK that is just the identity matrix times the number 
of stokeslets n: 

(K ) tJ = nSij (17) 

For the coupling tensor C, we have 

n 

Q=l 

n 

Co = ^rx], (18) 

a=l 

which is completely antisymmetric. If the origin is at the mean stokeslet position r* = 
iVP, then C = 0. 

Finally, when the body is rotating without translating, 

n n n 

(n 7)^ = r = XV x ( 7 v a )= 7 ^r"x (cu x f°) = 7 ^ ((r a ) 2 - f^r" • )cu, 

a=l a=l a=l 

so (Do is an inertia tensor: 

n 

{^h = Y,^ a ) 2 ^-< r D- ( 19 ) 

a=l 

12 



To second order, hydrodynamic interactions become important: 



v = v e + L(7v e ). 
Using our expression for the Oseen tensor gives 

87rr? ^ U f^)3 1 " i } 

(c ) 1 y r l***1* (^*^ ^ 



- ^ L ^ + 



87T77 ^ L r a/3 ( r a/3)3 

Since there is no term in the propulsion matrix which is zeroth order in 7, the expansion 
for T has the form 

T = T07" 1 + Ti + T 2 7 + • • • 

and will diverge as the effective stokeslet size approaches zero. Fortunately, the eigenvalues 
of T will not end up diverging as well. 

To see this, we will compute Y\_ by inverting £ in block form. We can identify the screw 
matrix S as the inverse of the Schur's complement of IK, giving S = (fi — C[K _1 C T ) _1 , and 
then T = — SCIK -1 . Likewise, the alacrity matrix A is the inverse of the Schur's complement 

of a 

When we neglect internal hydrodynamic interactions, the twist matrix then becomes 

T = -(fi " C Ko 'Cj^CoKo (20) 

As discussed above, there can be no twisting due to an applied force unless there are hy- 
drodynamic interactions. This means that the centers of twisting and reaction are the same 
here, and that at this point, Co = T = 0. As noted above, this point is also the mean 
stokeslet position r°. Since the twist matrix vanishes here, we can see that T will always 
have a null vector, regardless of where the origin is: to move the origin from the center of 
reaction to a position R corresponds to changing the twist matrix to J' = — So[-Rx]. But 
then J'qR = — So[-Rx]i? = 0, so J' still has at least one eigenvector, R, with a corresponding 
eigenvalue of zero. Thus as 7 decreases, one eigenvalue of T 7 _1 remains zero, though the 
other two may become large and complex. 
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When small hydrodynamic effects are added, the twist matrix expands to first order as 

T = Terr 1 + Ti. 

Because 7 is small, Ti makes a negligible correction to the To7 _1 term, except in the null 
space of To- Here, Ti? = J\R, which is independent of 7. To this order, the axis of spin is 

— * 

then R, the vector from the average stokeslet position r* to the forcing point. 

Since some eigenvalues of T can be complex, we cannot diagonalize it using real eigen- 
vectors. However, we can put T into Jordan canonical form using a basis of the form {v\, 
V2, R}- If we let i?J denote the dual of R, satisfying R d R = 1 and R d Vi = 0, then a real 
eigenvalue of T to this order will be 

A = R T d {l a /l + = R T d liR 

which is independent of 7. 

As noted above, the chiral response depends on hydrodynamic interactions between parts 
of the object. These interactions go to zero with the drag coefficient 7. Thus it is natural 
to anticipate that the angular velocity of the object should vanish with 7. Remarkably, this 
is not the case: we have just seen that a real eigenvalue of the twist matrix, and thus the 
angular velocity, reaches a non-zero limit as 7 — > 0. In this sense, there is a qualitative 
difference between the nearly free draining state and the perfectly free draining state. The 
difference may be understood through the propulsion matrix, which gives the force and 
torque in terms of the velocity V and angular velocity cu, and is regular as 7 — > 0. Both the 

— * — * 

amount of torque for a given V and no cu and the amount of torque for a given cu with no V 
are proportional to 7. With sedimentation, there is no net torque on the object, so we can 

— * 

find our cu for a given V by the requirement that the torque vanishes. If 7 is then reduced, 
both sources of torque are reduced in proportion, and the total torque remains zero with no 
change in cu. Thus cu has no tendency to vanish with 7. 

B. The tumble zone 

In Section III, we showed that the tumble zone had finite volume. The size and shape 
of this volume depend on the drag coefficient 7. We will now show that the volume of the 
tumble zone goes to zero at least as fast as 7 3 . Thus for sufficiently small 7, any collection of 
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stokeslets taken about any origin with no special symmetries will fall outside of the tumble 
zone, and must thus have simple fixed point chiral sedimentation. 

We use an argument similar to that in Section III, but choose the forcing point to be of 
the form R p = 'jQRp, where Q is independent of 7. About the center of twisting, we can 
write 

T = Ti 

S = So7 _1 + Si 

since in the low 7 limit, To = 0. Then 

T = Ti — (S0/7 + SOp^x] = Ti - Q(S + Si 7 )[HpX] (21) 

which has a part of order 7 and a correction of order 7 1 . The resulting discriminant A for 
the characteristic polynomial of T can be computed from Equation 12 as in Section III, but 
replacing S[i? p x] with QSoJ-RpX] +0(7(5). Using this substitution, we obtain a discriminant 
similar to Equation 13: 

1 r . n 3 



A = - 
2 



Tr[(QSo[i? p x])^] +om+0(Q b i) (22) 



Letting s 0m be the smallest of the eigenvalues of S gives the bound 

A > 2 S LQV + O(QS) + 0(Q 5 ), (23) 
except in the unphysical case that the stokeslets are perfectly collinear. In this case, one of 

— * 

the eigenvalues of So is zero, and taking R p perpendicular to this direction will make the 
Q 6 7° term vanish. 

Since 7 is small, the main contribution to the coefficient of the Q 6 term is from the 7 
part, which from Equation 23 is positive. For sufficiently large Q, we can then be assured 
that A > 0, giving one real eigenvalue for T. 

Thus we see that in the nearly free draining limit, the tumble zone can be fit inside of a 
sphere whose radius is proportional to the drag coefficient 7. As 7 — > 0, the tumble zone then 
must become vanishingly small. Unless the sedimenting object has the special property that 
its forcing point is exactly at the center of twist, we will thus get only one real eigenvalue 
for the twist matrix. We then expect globally stable chiral motion as it sediments. 
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VI. CHIRALITY 



The globally stable motion expected for nearly free draining sedimenting objects lends 
itself naturally to defining a chirality. If we denote the real eigenvalue of T by A, then 
A = Ld/F for F the magnitude of the applied force and u the constant angular velocity. We 
can try to use this A as a measure of the chirality. Conveniently, A is independent of 7 for 
nearly free draining objects, so we only need to know the shape of the object and the forcing 
point, and are not obliged to worry about the precise stokeslet strength. 

Unfortunately, if we try to use this measure to look for a "most chiral" object, we will 
be sorely disappointed: for a fixed 7, A diverges as the stokeslets become collinear. In 
this rather unphysical case, the eigenvalue of fl corresponding to rotations about the line 
of stokeslets will become zero, making £ non-invertible, and our expression for T, which 
depends on f) -1 , diverge. 



A. The distant forcing point limit 

In order to characterize the divergence of A we may simplify the analysis by considering 
the limit where the forcing point is far away from the stokeslets. This is a convenient choice 
because as long as the distance R p from the center of reaction to the origin is large, A is 
actually independent of the precise value of R p . This is true for any object, and does not 
depend on the approximation of small 7 used in Section V. 

To prove this assertion, we will first assume that we know the twist matrix around the 
center of twist. This choice of origin is somewhat arbitrary - any point close to the stokeslets 

— * 

will do. Once we have this T, we will move the origin to the point R p , where according to 
Equation 9 the new twist matrix is given by T' = T — S[i? p xJ. 

One of the eigenvalues of Sl-RpX] is zero. Since Tr(SlR p xj) = and Tr[(SlR p xj) 2 } < 0, 
the two nonzero eigenvalues must be imaginary. 

Next we will choose the basis, not necessarily orthogonal, which puts S[-R p x] into Jordan 
canonical form. Here, 



S[i2pX] = 



( T 1 \ 

r 
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where r is a generalized eigenvalue proportional to the pulling distance R p . We will define 
the basis {|0), |1), |2)} by S[,R p xJ|0) = 0, S[^,x]|l> = r|l), and $lR p x}\2) = r|2) + |1). 
We will also form the dual basis {(0|, (1|, (2|}, which satisfies = 5ij. 

Our goal is to find the real eigenvalue of ¥' = T — S[-R p x]. Since R p is large, T serves as a 
small perturbation of the S[i? p x] matrix. The real eigenvalue A must then be a perturbation 
of the single real eigenvalue of S[-R p x], namely zero. We will express its corresponding 
eigenvector as \v) = |0) + ei|l) + e 2 |2), choosing to scale it so that the coefficient of |0) is I, 
and €i <C 1. With this expansion, 

T» = X\v) 

T(|0) + e 1 |l) + e 2 |2))-S[ J R p x](|0) + e 1 |l) + e 2 |2)) = A(|0) + ei|l) + e 2 |2)) 
T|0) + eiT|l> + e 2 T|2)-eir|l)-e 2 r|2)-e 2 |l) = A(|0) + ei|l) + e 2 |2» 

The CjT terms must be small by comparison with the e^r terms, so we can drop them. Now 
applying (0| to both sides gives 

A=(0|T|0), (24) 
which is independent of the distance R p . 



B. Shape dependence of the chiral response 

Here we determine how the chiral sedimentation coefficient A depends on the locations 
of the stokeslets in the nearly free draining limit, in the case of distant forcing point. Even 
though A is independent of the distance to the forcing point in this limit, it can still depend 
on the orientation of the object relative to the pulling direction. We thus distinguish the 
coordinates of the stokeslets parallel and perpendicular to this forcing direction, denoted 
as z. We first note that our system has no distinguished origin, so A can depend only on 
the distances between the stokeslets. Accordingly, we measure stokeslet positions relative 
to their center, 




In terms of this, we define parallel and transverse radii of gyration, given by 

n 
a=l 
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FIG. 1: The stokeslet configuration used to check the scaling of A with the size of the object. To the 
left is a perspective view from an arbitrary direction. The pulling direction is toward the bottom 
of the cube, in the —z direction. The center and right views show projections of the stokslets onto 
the xy and xz planes, respectively. 



We use four parameters to characterize the distribution of stokeslets. The overall size 
can be expressed in terms of the radii of gyration given above. In addition, we use a length 
Z denned below to characterize inhomogeneity in longitudinal position, and a dimensionless 
quantity A to characterize anisotropy in the transverse plane. 

To simplify matters, we will focus on configurations with the fewest number of stokeslets 
required to make a chiral response possible. Since the object as a whole also includes a 
forcing point, we only need three stokeslets to guarantee a non-planar configuration. In 
such cases, with a distant forcing point, there are nine coordinates which can specify shape. 
However, A is independent of translation and of rotation around the pulling axis, so only 
five coordinates are potentially significant. We next show that the four parameters named 
above appear to suffice. 

To begin, we check the dependence of A on the size of the object. We do this by fixing a 
configuration of stokeslets and then computing A as we uniformly change the inter-stokeslet 
distances. The particular configuration we use is shown in Figure 1. It has the three 
stokeslets arranged so that their projection in the xy plane is an equilateral triangle with 



and 




a=l 



The total radius of gyration is then R 2 g — R?, + R\ 
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(a) (b) (c) 




R R l R ll 



FIG. 2: Numerical results showing the scaling of A with R (a), R± (b), and i?y (c), for configurations 
like that of Figure 1. In (a) there is a single scaling exponent of -2. In (b), for R± 3> R\\ = 1, we 
have a scaling exponent of -5, and for R± <C R» = 1, it is constant. In (c), we have an exponent 
of 3 for #|| «fli = l, and -2 for i?y » R± = 1. 

side length i? centered about the origin, and their positions along the £ axis are and ±R. 
In this case, it does not matter which corner of the triangle is at which z value; by symmetry, 
rearranging them can at most change the sign of A, while its magnitude is our concern here. 
We will then move the forcing point to R p z, with R p ^> R, and compute A. 

In Section II we noted that the propulsion matrix depends linearly on 77, so A ~ rf 1 . 
We can ignore this simple dependence on viscosity by setting 77 = 1. We will also set 
7 = Q-K7]p = 1CT 2 , with p in the same arbitrary distance units we use to measure R. As long 
as p <C R, this is within the regime of small 7, so the precise value does not matter. 

As discussed in Section II, all elements of T scale as an inverse length squared, so A must 
as well. Since A is independent of the Stokes radius and distance to the forcing point, as 
shown in Sections VA and VIA, we must form this length scale from the inter-stokeslet 
distances. Indeed, we can verify numerically that A ~ R~ 2 , as shown in Figure 2(a). Since 
R g scales with R, it is clear that A ~ R^ 2 - We can further try to break this dependence 
down into one based on R± and R» . To begin, we fix the z positions of the stokeslets to be 
and ±1, and then vary the side length of the equilateral triangle. As shown in Figure 2(b), 
when the side length is long compared to the z positions, we get A ~ R~]_ 5 . This corresponds 
to a flat transverse object. For small side lengths, we get A ~ R ± . This corresponds to an 
object that is elongated along the pulling direction. We can also see what happens when we 
fix the side length of the equilateral triangle in the transverse projection at 1, and instead 
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vary the z distance between stokeslets, putting them at and ±R\\. The results are shown 
in Figure 2(c). We see that for i?y <C 1, we get A ~ and for i?y 3> 1, we get A ~ i?^ 2 . 
Taken together, these observations suggest that we can write A = R~ 2 f(Ri\/R±), where 



m~< : (25) 




We can see that when R\\ 3> R±, the function / is a constant. Thus in this regime we know 
the scaling of A based on relative transverse and longitudinal sizes, and can focus on other 
aspects of the object's shape. 

We will consider two general distortions of our shape from the previous one: first, we 
will relax the requirement that the z values be equally spaced, in order to see the effect 
of bunching a pair of stokeslets together. Next, we will remove any restrictions on the 
transverse shape. 

To characterize the bunching, we will use the inverse squared moment Z, defined by 



n 

c\-2 



1 - r z 

a=l 



This length Z is dominated by the closest pairs of stokeslets. If we consider the ratio Z/R\\, 
we get a dimensionless quantity which becomes large if some stokeslets are bunched close 
together. 

If a pair of stokeslets is bunched together, the hydrodynamic interactions between them 
become stronger. We expect this greater interaction to promote chiral behavior. Indeed, 
our numerical studies indicate that uneven spacing leads to larger A. We again fix the 
transverse projection of the stokeslets to be an equilateral triangle, with side length a. We 
then choose the two extremal longitudinal projections to be at ±5, and allow the middle 
stokeslet position z to vary between the other two. Figure 3 shows the chirality as a function 
of z, for three different values of a. We see that there is a peak in A as the stokeslets approach 
each other, but it falls off if they get too close. The maximal A occurs when the longitudinal 
spacing is about equal to the transverse spacing. 

To characterize the shape of the transverse projections, we consider the eccentricity of 
the inertia ellipse. If we define a projected tensor of inertia by 



n 

a=l 



1 U 

-£(r?-r?)(r?-r?) 
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FIG. 3: The chiral coefficient A for a three-stokeslet object whose transverse projection is an 
equilateral triangle of side length a. Two of the longitudinal coordinates are fixed at ±5, and the 
third is varied over z values between them. 



for i,j e {x,y}, then we can use 



4Det(l) 
Tr 2 (1) 



as a measure of the eccentricity. It goes to zero when the stokeslets are collinear, and one 
when they are isotropically arranged. 

We can now consider A as a function of both A and Z/R\\. We confine ourselves to shapes 
with R p ^> R\\ ^> R±, which gives maximal A as seen above. 

To see the dependencies, we generated 10 4 random 3-stokeslet configurations, choosing 
each stokeslet from the box [—1/2, 1/2] x [—1/2, 1/2] x [—10, 10], and discarding it if R\\ < 
10R±. Again we removed 7 and 77 dependencies by taking 7 = 10~ 2 and 77 = 1. 

The observed A values varied widely and irregularly. However, if we define the pth moment 
of the stokeslet positions 

n n -1 1/p 

EE 



u p = 



n(n — 1) 



If** — fP\p 



a=l @=a+l 

and instead plot A(f/_ 2 ) 2 (Z/_R||) 2 , we get a relatively smooth bounded function. Thus we 



can write 



A 



gun) 



(26) 
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FIG. 4: The g function plotted versus A and {Z/R\\) 
high A and Z/R\\ . 



2 . We can see that it is bounded and prefers 



where g is a bounded function of its arguments. 

Figure 4 shows a plot of g as a function of A and (Z/R\\) 2 . From this plot, we can see 
a definite dependence on A, indicating that g, and thus A, prefer higher A. This means 
that faster rotation occurs when the transverse projection is isotropic rather than elongated, 
while the object as a whole is long and slender. 

In general, studying this simple 3-stokeslet case in the limit of distant forcing points has 
shown that the preferred shape for high chirality is a long and slender object. Along the 
length of the object, some clustering of stokeslets is preferred, and in the transverse plane 
it is beneficial to have an isotropic arrangements of stokeslets. 

So far we have only considered the magnitude of the chiral response for our three stokeslet 
systems. It would be convenient if there was an easy way to determine the sign of the chirality 
as well. We propose a method which seems to give acceptable results for those systems with 
large values of |A|. 

We first order the stokeslets according to their longitudinal proximity to the forcing point. 
In the transverse projection, the ordering will form either a clockwise or counterclockwise 
triangle. We propose that these respectively correspond to a negative and a positive chirality. 
The physical argument for this triangle rule is that as the object sinks, the first stokeslet 
will have a stronger interaction with the second than the third, and so on. This will cause a 
slipstreaming effect, where the fluid behind the first causes less drag on the second behind 
it. This preferentially allows the object to move in that direction, much like a corkscrew. 

To test this numerically, we generated 10 4 triples of stokeslets chosen at random from 
the box [—2,2] x [—2,2] x [—2,2]. For each object we computed the chirality and applied 
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All configurations 


Very chiral configurations 


Slender configurations 


Matching signs 


7518 


2498 


9868 


Different signs 


2482 


264 


162 


Percent matched 


75.2 


90.4 


98.7 



TABLE I: Comparison between the number of times the actual sign of the chirality matched the 
sign estimated using our triangle rule. The very chiral configurations were selected from the rest 
via the criterion that their chirality be larger than the average. In addition, 10 4 configurations 
were chosen at random from the slender regime studied earlier, where we expect to find the most 
chiral configurations. 

the above triangle rule. The results are shown in Table I. The triangle rule predicted the 
correct chirality roughly three quarters of the time. We anticipate as well that more chiral 
objects will be more likely to follow our sign convention, as the slipstreaming effect will 
be stronger. To test this, we repeat our comparison using only those configurations whose 
|A| value was larger than the average. As shown in Table I, our method was indeed more 
accurate with the more chiral configurations. We can test this in another way by limiting 
ourselves to the more chiral configurations which we know arise when our three stokeslets 
are instead chosen from the box [—1/2, 1/2] x [—1/2, 1/2] x [—10, 10]. In the case of these 
slender configurations, our method is quite effective. While it is not perfect, it can provide 
a reasonable guess at the sign. 

Our explicit calculations above focused on the simplest stokeslet object that can have 
chirality: three stokeslets with the forcing point at infinity. We noted that such an object 
has five relevant degrees of freedom, but studied the effect of only four of them. To specify 
the minimal object completely therefore requires an additional parameter. One choice is to 
use the full 3x3 inertial tensor, instead of its transverse projection. The principal axes of this 
tensor need not be aligned with the forcing direction, so we can take our additional parameter 
to be the smallest angle between a principal axis and the forcing direction. Evidently for 
the elongated objects with large A we have been studying, this angle is small, and does not 
have a major effect in this regime. 
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VII. EXAMPLES OF BEHAVIOR 



A. Numerical results 

As a simple test of the results from Section V, we can generate several stokeslet configu- 
rations at random, and verify that in the nearly free draining limit we get the simple chiral 
sedimentation predicted above, with the expected axis of rotation and angular velocity. We 
will do this with four objects: For object A, we form a five stokeslet object by picking 
random positions in the box [—2, 2] 3 and setting the origin as the forcing point. Object B 
is the same as object A except for the location of the forcing point. This point is moved 
closer to the center of twisting in order to increase the tendency to tumble. Specifically, 
the center of twisting is determined at a particular choice of stokeslet radius, namely 2/3 of 
the radius pmax which would create contact between stokeslet spheres. The forcing point 
is then placed at this center of twisting, and remains there as the stokeslet radius is varied 
and the resulting motion measured. Object C is created the same was as object B, but with 
a different random choice of stokeslet positions. Finally, object D is a random ten stokeslet 
object, again with the origin moved as above. These are shown in Figure 5. 

For each object, we first determine the axis we expect the object to rotate around in the 

— * 

nearly free draining limit. This is easy: as described in Section VA, the real eigenvector A 
of the twist matrix is just the vector from the forcing point to the average stokeslet position, 
X = r°. 

To find the angular velocity, we compute T and Ti as in Section V, and form the basis 
vi, V2, and A which puts T into Jordan canonical form. Let A J be the dual of A, which 
satisfies AjA = 1 and AjVj = 0. Then the real eigenvalue of T is just A = AjliA, and we can 
find the angular velocity from uj = X\F\ with \F\ the magnitude of the sedimenting force. 

In Figure 6, we compare these nearly free draining results with the results obtained from 
inverting Equation 16, using \F\ — r] — 1. We see that the nearly free draining results hold 
over several decades of 7 values. Significant deviations occur only when the object is near 
the tumble zone. In the tumble zone, there is no single value of uj or cos 6 which can be 
plotted. However, we can see that we need to be quite close to the center of twisting for 
this to occur; the global stability and predictions from the nearly free draining limit made 
earlier are quite robust in practice. 
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►X 

FIG. 5: Projections of the four stokeslet configurations used in Section VII A. In each image, the 
grey circles represent stokeslets, a small square marks the average stokelet position, and arrows 
point to the forcing points used. Objects A and B are identical except for the positions of their 
forcing points, and are shown on the left. Objects C and D are shown in the middle and on the 
right. Each object has been rotated so that the coordinate axes are aligned with the principal axes 
of the inertia tensor, with z and x corresponding to the largest and smallest of these, respectively. 
The size of the grey circles corresponds to the largest the stokeslets can be without causing the 
object to enter the tumble zone. In the case of the leftmost images, this is done with respect to 
object B. 

We next study the effect of initial orientation on the sedimenting path, as Makino and Doi 
did for their skew propeller shape [3]. We do this by taking N = 100 objects, each shaped 
as the object plotted in Figure 6(a) above, but with different random initial orientations. 
We then release them from the same point (x,y,z) = (0,0,0), and consider their positions 
as functions of time, ignoring interactions between different objects. We determine these 

— * — * — * 

positions from the velocities given by V = A(t)F and cu = J(t)F. A(0) and 1(0) can be found 
from Equation 16, and their time evolution is governed by the differential equations given 
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FIG. 6: A comparison between the predicted values of the axis of rotation and angular velocity in 
the nearly free draining limit with the full results valid for all values of the drag coefficient. On 
the left vertical axis, the solid line represents the full u> when the object is undergoing the globally 
stable chiral motion, and the dotted line represents the value computed from the perturbative 
expressions in Section V. On the right vertical axis is the cosine of the angle 6 between the axis 
of rotation and the axis of rotation computed in the nearly free draining limit, (a), (b), (c), and 
(d) correspond to objects A, B, C, and D, respectively. In the case of (a), the object does not 
enter the tumble zone at all; before this happens, p has increased to the unphysical point where 
the stokeslets overlap. However, the rest of the objects had their forcing points chosen in a manner 
which required them to be in the tumble zone for larger values of p, and their plots break off before 
p max is reached. 



in Section II. While the elements of the these matrix are coupled together, the equations 
have no singularities, so we use Mathematica's NDSolve function [14] to numerically find 
the solutions and expect reasonable accuracy. We use rj = 1, 7 = 10~ 2 , and supply a force 
F = z. 
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FIG. 7: The skew propeller shape used by [3]. The two orthogonal disks have radius a and are 
fixed so their centers are a distance 2£ apart. The center of twisting for this object coincides with 
the origin, so its twist matrix has three real eigenvalues. 



These results can be compared to those from the skew propeller shape, as well as a simple 
ellipsoid. The propeller consists of two orthogonal disks of radius a attached via a thin rod 
so that their centers are a distance 2£ apart, as shown in Figure 7. The relevant portions of 
the mobility matrix are A =diag(a :r , a x ,a z ) and 



T 



^0 60^ 
6 




with 



Ox = 



3 (4a 2 + 5£ 2 ) 
64ar/(5a 2 + U 2 ) 
3 



64ar) 



3£ 



64ar/(5a 2 + 6£ 2 ) 

In the following, we use £ = 3a, and then set a = 1 in the same length units we used for our 
nearly free draining object. 

The skew propeller is an example of an object whose twist matrix is symmetric, and thus 
allows us to compare our nearly free draining object with something in the tumble zone. 
The ellipsoid allows a comparison with an object that has no translation - rotation coupling; 
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FIG. 8: Plots of normalized w(t) and h(t) as functions of time. The width increases linearly 
with time for the ellipsoidal particle, but remains bounded for the particles with a nonzero twist 
matrix. The spread of the particles increases linearly with time for both the skew propellers and 
the ellipsoids, but after an initial transient remains constant for the nearly free draining particles. 

its twist matrix is zero. We will choose its dimensions so that its alacrity matrix is the same 
as that of the skew propeller. 

To do the comparison, we can look at the width of the distribution of particles as a 
function of time, as well as the spread in the z direction: 

1 N i 

w ^ = jvSv^ 2(t) + ^ (t) (27) 

i=l 
1 N 

h(t) = -<*(*)>! (28) 

i=l 

where {xi(t), Ui(t), Zi(t)} is the position of the particle at time t, and (z(t)) = jj Yli z i(f) ls 
the average z position of the ensemble at time t. 

Figure 8 shows w and h, normalized by the maximum linear distance between two points 
on the object, l m . The ellipsoids must distribute themselves on the surface of a sphere 
sinking at a constant velocity whose radius increases with constant velocity [3] . Thus h and 
w are both linear in time for ellipsoids. 

The widths of the distributions for the skew propellers and our sample object are both 
bounded. The skew propellers evidently approach a constant w, while the nearly free drain- 
ing objects have a w which oscillates at their rotation frequency u. The spread h for the 
skew propellers in the tumble zone increases linearly. However, we see that after an initial 
transient motion, the longitudinal spread of our sample particles remains constant. This is 
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because all of them begin to sediment in the same regular manner. 

Thus, overall we see that the ellipsoids spread out into a spherical volume as they sink, 
with radius increasing linearly. The skew propellers spread out in a cylindrical shape parallel 
to the applied force. The length of the cylinder increases linearly, while the radius undergoes 
decaying oscillations about a value smaller than the linear extent of each object. The nearly 
free draining particles spread out over a flat disk with constant longitudinal spread and a 
radius which oscillates at the same frequency that each particle spins at. The amplitude of 
these oscillations is slightly larger than the maximum extent of the object. 

B. Experimental illustrations 

In order to verify that the chiral rotation discussed above is significant in practice, we 
created some arbitrarily shaped bodies and observed their sedimentation. We used both 
viscous and non-viscous solvents. This allows us to gauge the importance of inertial effects. 

For the viscous solvent, we cut small objects out of a rod of nylon plastic, a few millimeters 
in length in their longest direction. We also took small lengths of copper wire and bent them 
into twists or knots. Our objects were allowed to sediment in a 700 mL beaker filled with 
vegetable oil. Such oils have kinematic viscosities of the order 30 cSt [15], and our nylon 
pieces fell at around 0.2 cm/s, giving a Reynolds number of slightly less than 1, well within 
the Stokes regime. The copper twists, which fell more quickly, are still at Reynolds numbers 
where inertial effects are not expected to be important. 

We used tweezers to hold each object just below the surface, then released it and used a 
camera to take pictures at a rate of about 3 frames per second. For these uniform materials, 
the forcing point is the center of mass, which we expect to be close to the center of reaction. 
Thus it is not clear from our arguments above that these objects should be outside the 
tumble zone. Nevertheless, we were able to see chiral sedimentation with many of these 
objects. Figures 9(a,b) show multiple-exposure views for both a nylon piece and a twist of 
fine copper wire. The helical path is obvious for the copper piece, but less so for the nylon. 
Figure 9(c) shows the same nylon piece, in a separate run, from above. Here the helical 
nature of the path is easier to see. 

In addition to the objects shown, we tested over a dozen other objects made in the same 
way. Some displayed little or no rotation, and simply settled into a preferred orientation. 
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FIG. 9: (a): A multiple-exposure image of an irregular piece of nylon sedimenting in vegetable oil. 
The nylon piece is a few millimeters in length, and the pictures were taken about 1 second apart. 
It is clearly rotating around the vertical axis, (b): A multiple-exposure image of a fine piece of 
copper wire sedimenting in vegetable oil. The pictures were taken about 0.3 seconds apart, and 
the distance scale is the same as in (a). The object is rotating about the vertical axis as it follows 
a helical path down, (c): A multiple-exposure image of the same piece of nylon in (a), though not 
at the same time. From above, the helical path is more apparent. More images and movies are 
available at http://jfi.uchicago.edu/~tten/ChiralSedimentation/ 

Some of the heavier ones, which sank very rapidly, showed a slight rotation about axes other 
than the vertical. We cannot tell if this was actually a case of the objects tumbling; we 
suspect that it was instead an initial reorientation which aborted when they hit the bottom 
of the beaker before reaching their preferred orientation. 
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For comparison, we dropped small shards of brittle plastic into a salt water solution, 
whose viscosity was lower than the oil's. The objects, cut from a disposable spoon, were 
a few millimeters in size. Salt was added to the water to achieve nearly neutral buoyancy 
without greatly affecting the viscosity. In this solution, the objects fell at around 1 cm/s, 
giving a Reynolds number < 100, which is not fully in the regime of Stokes flows. However, 
we still observed chiral sedimentation, so even at this Reynolds number the inertial effects 
do not appear to change the motion qualitatively. 

In these studies, we monitored for residual circulation in the water by putting a small 
cylinder of floating plastic on the surface. This cylinder remained stationary, indicating that 
any residual flow is much smaller than the chiral motions. 

With these plastic pieces, no ongoing tumbling motion was seen; either there was no 
rotation, or else they rapidly reoriented themselves and twisted around the vertical axis. 
Figure 10(a) shows a multiple-exposure picture of a typical path. This object, shown close- 
up in Figure 10(b), turned to the same preferred direction regardless of initial orientation, 
and always rotated with the same sign. However, this is not the only behavior; the object 
pictured in Figure 10(c) had two opposite orientations which were stable. These produced 
opposite signs for the rotations. In addition, some objects showed negligible rotation, though 
they did go to the same stable orientation. This must correspond to an instance where A ~ 0. 

VIII. DISCUSSION 

In the foregoing we have explored how slowly-sedimenting noncompact objects of generic 
shape rotate as they sink, revealing chiral structure. These objects were represented as 
collections of stokeslets, which are known to provide a good representation of a broad range 
of real objects [11]. We infered the propulsion matrix from the matrix of Oseen interactions 
between pairs of stokeslets. This propulsion matrix is sufficient to determine the entire 
motion under slow sedimentation at low Reynolds numbers [9]. To determine the chiral 
rotation, it is sufficient to know the 3x3 twist matrix T derivable from the propulsion 
matrix. In the case when T has only one real eigenvalue, there is globally stable motion 
corresponding to rotation about the corresponding eigenvector [5]. 

Though all chiral rotation must vanish when there are no hydrodynamic interactions, 
in the nearly free draining limit where these interactions are arbitrarily small, there is 
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FIG. 10: (a): Three multiple-exposure pictures of the 9 millimeter long object pictured in (b) as it 
sediments in salt water. The pictures were taken about 0.15 seconds apart. In each case the twist 
about the vertical axis as it moves in a helix is clearly visible, indicating that the chiral effects on 
sedimentation are similar to those in the viscous solvent of Figure 9. Each picture corresponds to 
a different initial orientation of the object. Though the transient motion was different in each case, 
it always ended up turning to the same preferred orientation and twisting in the same direction, 
(c): Another object cut from a plastic spoon. This object has two stable orientations, which lead 
to twists about the vertical axis in opposite directions. 

nevertheless a constant and finite rotation about a fixed axis. The angular velocity in this 
limit is independent of the strength of the interactions, and the rotation axis approaches the 
line between the forcing point and the center of reaction. 

The features of an object that determine its chiral sedimentation are unexpectedly sub- 
tle. Indeed, the rotation rate depends on the stokeslet positions in a singular way, with 
unevenly spaced stokeslets giving the largest response. For such configurations it is the 
nearest distance that dominates, and small displacements of the stokeslets on the order of 
this shortest distance suffice to reverse the sign of A. Thus A is not a gross indicator of 
overall chiral shape. Instead, it is a local probe, sensitive to local orientations relative to 
the overall object. The maximum responses occurred for thin, screwlike objects. Similar 
objects at the microscopic scale include biological filaments such as f-actin or microtubules. 

The connection between our simple stokeslet objects and real objects has not been fully 
explored in this paper. Carrasco and de la Torre [11], for example, describe methods for 
implementing the stokeslet model which appear applicable to the objects we discuss. Thus, 
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rather than predicting the chiral response of any real object, we focused instead on finding 
the scaling and analytical asymptotic behavior for nearly free draining objects. 

We have developed an empirical rule to predict the sign of the chirality for some simple 
objects. However, this method should be improved. We would like to find a simple method 
to determine the chiral sign that is not only more accurate, but will also generalize to 
arbitrary objects. We also would like to establish analytically the scaling that we empirically 
determined in Section VI B, and to include the effects of brownian motion. 

The free draining limit we use is physically approachable for the sedimentation of some 
large molecules or other polymers, formed by assembling macromolecules or colloidal parti- 
cles. One could conceive of attaching a fluorescing group to such a molecule, and then using 
fluorescence polarization in a centrifuge to measure the spinning rate. The spinning rate 
could be used to characterize the object. 

Even in cases where the objects are not nearly free draining, we expect most of our 
conclusions to apply qualitatively; the nearly free draining limit is not the only way to 
escape the tumble zone, and general objects without symmetry will often see the globally 
stable behavior. 

The chiral sedimentation treated here is only one example of how a colloidal object of 
irregular shape might respond in a chiral way. For example, objects sedimenting in shear 
flows can undergo net lateral drift according to chirality [16]. Varying sedimenting forces 
periodically in time could also be used to probe further properties of the propulsion matrix. 
Molecules of sub-micron scale such as folded RNA must also exhibit chiral sedimentation, 
though they will be greatly influenced by thermal brownian motion. Beyond the context 
of hydrodynamics, such objects can show chirality via their self-assembly properties. For 
example, two copies of a chiral globular protein have a most favorable orientation for binding. 
When many such copies self-assemble in this way, the least constraining mode of assembly 
is a one dimensional stack. Such a stack must in general show a chiral twist which may 
limit the stack's potential to stick to its neighbors. Aggeli et. al. [17] use this as a model 
for the formation of peptide fibrils. This generic view may account for the prevalence of 
one-dimensional assemblies of biological molecules. Such responses are a promising course 
of study for the future. 
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IX. CONCLUSION 



The most classic chiral response of microscopic matter, the rotation of the polarization of 
light, has been studied for over a century. Here we have discussed an equally fundamental 
response: the chiral interaction of an irregular object with a surrounding viscous liquid. 
In this case the chiral properties arise entirely from the object's geometry. We have seen 
that macroscopic objects of arbitrary shape have readily observed chiral sedimentation. The 
greatest response seems to occur when the drag is concentrated at one end of an elongated 
object. This study is only a first step towards understanding how shape creates chiral 
responses in colloid-scale materials. There are numerous ways to explore various shapes and 
numerous other responses, as sketched above. Understanding how shape determines chiral 
response should be valuable as a way of assessing the shapes of unknown objects and as a 
way of designing shapes to create desired responses. 
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